A PROOF FOR THE INFORMATIONAL COMPLETENESS OF THE 
ROTATED QUADRATURE OBSERVABLES 



J. KIUKAS, P. LAHTI, AND J.-P. PELLONPAA 

Abstract. We give a new mathematically rigorous proof for the fact that, when S is a dense 
subset of [0, 2tt), the rotated quadrature operators Qg, 6 S, of a single mode electromagnetic 
field constitute an informationally complete set of observables. 



1. Introduction 



Since the pioneering works of Vogel and Risken [15] and Smithey et al [13] the measurement 
of the rotated quadratures Qg, 9 G [0, 27r), has formed one of the major tools in the quantum 
state tomography, see, e.g., the compilation [TT]. The basic idea behind the state reconstruction 
is well known: the inverse Radon tranform of the quadrature measurement statistics allows one 
to reconstruct the Wigner function of the state in question and the Wigner function separates 
states, so that the statistics uniquely determines the state. We do not question the validity 
of this argument. However, the existing literature, which we are aware of, does not give a 
full justification that this procedure actually applies to all possible states of a quantum system. 
Therefore, in this paper, we wish to give a direct proof of the fact that the set of quadrature ob- 
servables Qe, G S, S C [0, 27r) is dense, is informationally complete, that is, the measurement 
outcome statistics p® e of these observables uniquely determine the state T. 

There is a beautiful group theoretical proof of the informational completeness of the observ- 
ables Q e , 9 G [0, 27r) [5]. This proof builds on a general method of constructing informationally 
complete sets of observables using the theory of square-integrable representations of unimodular 
Lie groups. The result in question is then obtained as a special application of this theory to the 
Weyl-Heisenberg group. Due to the practical importance of the result, we give an alternative 
direct proof of it. The proof of this fact in Section [3] forms the main body of this paper. 

In the final Section [4] we comment on the measurability of the quadrature observables and 
we compare the tomography based on their measurements on the one obtained from the phase 
space measurements. 



2. Basic notations and definitions 

Let TC be a complex separable Hilbert space, and L(7i) the set of bounded operators on TC, 
and T(TC) the set of trace class operators. We let || ■ ||i denote the trace norm of T(TC). (The 
operator norm of L(TC) will be denoted by || ■ ||.) When TC is associated with a quantum system 
(such as the single mode electromagnetic field), the states of the system are being represented 
by positive operators T G T(TC) with unit trace, density operators, and the observables are 
associated with the normalized positive operator measures defined on the Borel a-algebra B{R) 
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of the real lineQ Among them are the conventional von Neumann type of observables, that is, 
projection valued measures P : B(JRL) — > L(TC), or, equivalently, selfadjoint operators in Ti.. 

The measurement outcome statistics of an observable E : B(M) — > L(H) in a state T is given 
by the probability measure X i— > ti[TE(X)] —: pf (X). 

Definition 1. A set M. of observables E : £>(R) — > L(H) is informationally complete, if any 
two states S and T are equal whenever tr[TE(X)] = tr[SE(X)] for all E e M and X G £>(R). 

Thus, the informational completeness of a set M. of observables means that the totality of 
the measurement outcome distributions p?, 6 M, determines the state T of the system. 
Clearly, a set M. of observables is informationally complete if and only if T = whenever T is 
a selfadjoint trace class operator with tr[TE(X)] = for all E G M. and X G B(R). We will 
use this characterization in our proof. 

Fix {|n) | n G N} to be an orthonormal basis of H. (This is identified with the photon 
number basis in the case where Tt is associated with the single mode electromagnetic field.) 
Here N = {0, 1, 2, . . .}. We will, without explicit indication, use the coordinate representation, 
in which TC is represented as L 2 (R) via the unitary map Ti 3 \n) i— ► h n G £ 2 (R), where h n is 
the nth Hermite function, 

1 _1 2 

h n (x) = = H n (x)e z x , 

and H n is the nth Hermite polynomial!! Let a and a* denote the usual raising and lowering 
operators associated with the above basis of H; they are consider as being defined on their 
maximal domain 



D(a) = D(a*) = |y? G Tt\ ^Tn\(<p\n)\ 2 < oo 



n=0 



Then define the operators Q = -^(a* + a) and P = ^(a* — a), which, in the coordinate rep- 
resentation are the usual multiplication and differentiation operators, respectively: (Qip)(x) = 



xip(x), (Pip)(x) = —i-£(x). (Here the bar stands for the closure of an operator, so that e.g. Q 
is the unique selfadjoint extension of the essentially selfadjoint symmetric operator -^(a* + «); 
see [12, Chapter IV] or [H Chapter 12] for details concerning the domains of these extensively 
studied operators.) In the case of the electromagnetic field, Q and P are called the quadrature 
amplitude operators of the field. In addition, a = ^{Q + iP) and a* = -^(Q — iP) (see 
e.g. [121, p. 73]). The Schwartz space <S(K) of rapidly decreasing C°°-functions is included in 
D{Q)r\D{P) = D{a). 

For a function g : R — > R, we let g( k > denote the kth derivative of g (with g(°> = g), provided 
it exists. We need the following elementary commutation relations, which hold whenever if G 
iS(R), and g : R — ► R is continuously differentiate and bounded: 

(1) (g(Q)a*-a*g(Q))<p = ^ {1) (QV, 

(2) (g(Q)a - ag(Q))<P = -j= 2 9 {1) (Q)v- 



formalized positive operator measure is a map E : i3(R) — > L(H) which is er-additive in the weak operator 
topology, and has the property E(W) — I, that is, for which X i— » tr[TE(X)} is a probability measure for each 
positive trace one operator T . 

2 Hermite polynomials are, of course, given by the following recursion relation: Hq{x) — 1, H\(x) — 2x and 
H n+ i{x) = 2xH n (x) - 2nH n -i(x). 
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Let N denote the operator a*a. It is selfadjoint on its natural domain 

oo 

d(n) = (tp g n\ J2 n2 M n )\ 2 < °°}- 

n=0 

Now the phase shifting unitary operators are e 10 , and we can define the rotated quadrature 
observables Qe by 

Qe = e l6N Qe- i6N , 9 G [0, 2ir). 
The spectral measure of Qe will be denoted by P Qe : <B(R) — > L(H). 

3. The proof 
We need the so called Dawson's integral 



daw(x) = e~ x ' 2 / e* 2 dt 
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(see e.g. [U pp. 298-299] or [HI Chapter 42]). The following lemma lists those properties of 
Dawson's integral that we are going to use. Since the Dawson's integral has been studied 
extensively, they are probably well known. However, as we were unable to find these results 
directly stated in the literature, we give a proof in the Appendix; a reader familiar with the 
results may skip that proof. 

Lemma 1. (a) daw : R — > R is a C°° -function, and 



GO 



lim daw (fc) (x) = for all keN. 

x— >±oo 



1 °° f 11^ ! 
daw (1) (x) = - , ' H 2n (x) for all x G 

rc=0 v 7 



where H n is the nth Hermite polynomial. 

Lemma 2. There exists a C°° -function f : R — > R, swc/i £/m£ 

(i) eac/i derivative p n > of f, n = 0, 1, 2 . . . is a bounded function, and 

(ii) (n|/(Q)|n) = 5„0; /or all n G N. where 5 is the Kronecker delta. 

Proof. We put / = 2daw ( - 1 \ According to Lemma [U / is a C°°-function and (i) holds. Lemma 
ED (b) gives the expansion 

(3) /(*)=zy^(*), *gr. 

Since the series in ([3]) converges pointwise, we get 

00 (—'\) n ri \ 00 (—i^l 



for each x G R. In addition, ^^^(nl) 2 / (2n)\ < oo, so the series converges also in L 2 -norm. 
Noting also that H 2 h G £ 2 (M), we can justify the following computation. 

(hn\f(Q)h n ) = —{HZhotfho) = —^^^{H^holV^m^) 



k=0 K ' 



1 ^ (-l) k k\ 2 k+n y^(2k)\(n\) 2 _ ^ (-l) fe n! 
2^! ^ 2 fc (2fc)! v /* : (n - k)\{k\) 2 ~ ^ (ra- fc)!fc! 



1 — n \ / 



J n0- 

\ K I U-*i 

k=0 

Here the finite sum after the third equality is obtained by noting that h 2 k is orthogonal to 
H 2 h whenever k > n, which is due to the fact that the latter function is a linear combination 
of Hermite functions ho, ... , h 2n . The fifth equality follows from the formula 7.375(2) of p. 
838]. □ 

Now we choose a function / satisfying the conditions of LemmaEl The function / will remain 
fixed throughout the rest of the paper. 



Lemma 3. (n + k\f^{Q)\n) = (-l) k V¥k\5 Qn for all k,n G N. 

Proof. We proceed by induction with respect to k; the initial step is provided by condition (ii) 
of Lemma El The induction assumption is that for some k G N, the equality 



(n + k\f k \Q)\n) = {-l) k V¥k\5 Qn 
holds for all nGN. We must show that 

(n + k + l\f^ k+1 \Q)\n) = (-l) k+ W^ k+1 (k + l)\5 0n , n G N. 
But by using (pQ), we get 

(n + k + l|/( fc+1 )(Q)|n> = y/2(n + k + l\f k \Q)o*\n) - y/2(n + k + l\a* f< k \Q)\n) 
= y/2(n+l)((n + 1) + k\f k \Q)\n + 1) - y/2(n + k+l){n + k\f^ k \Q)\n) 



= -^2{n + k+ l)(-l) k V¥k\5 an = (-l) k+ W^ +1 (k + l)\5 0n 

for all n G N by the induction assumption. (Note, in particular, that the first term in the 
expression following the second equality is indeed zero by the induction assumption, because 
71 + 1 >0.) □ 

Lemma 4. (n + jfe|/(*+»)(Q)|n) = 2 l (-l) k y/2H\(l + k)\ 5 ln , k,l,neN, n>l. 

Proof. Now we use induction with respect to I, so the initial step / = is given by Lemma El 
The induction assumption is that 

( n + k\f k+2l \Q)\n) = 2 l (-l) V2 fc /!(/ + k)\ 5 ln 



holds for some / G N, all k G N, and all n> I. We have to show that this holds also when / is 
replaced by I + 1. Accordingly, let k G N, and n G N with n > I + 1. Using the commutation 
relation (j2j) and the induction assumption, we get 

(n + k\f k+2 ^\Q)\n) = -V2(n + k\f k+2l+l \Q)a\n) + y/2(n + k\af (k+2l+1) (Q)\n) 
= -y/2^{{n -l) + (k + l)\f^ +2l \Q)\n - 1) + \j2(n + k + l)(n + (fc + l)|/« fe+1 ) +2 ')(Q)|n> 
= -V2^2\-l) k+1 ^2 k +H\(l + k + l)\ = (-1) V2(/ + 1)2^^2*^ + 1 + fc)! <$, +1>n 
= 2' +1 (-l)V2 fc (/ + !)!(/ + l + fc)!<W 

Here the induction assumption was applied to the first term following the second equality with 
n and k replaced by n — 1 and k + 1, and to the second term with n and replaced by n and 
k + 1. Here it is essential to note that n > Z + 1 > I, which makes the second term zero. □ 

In order to construct the proof for the informational completeness of the quadratures, we 
need some additional tools. First define, for each fixed fceN and X G <6(R), 

"27T 



1>2TT 

V k (X) := / e- ike P Qe (X) d6 G L{H), 
Jo 



where the integral is to be understood in the cx-weak operator topology. Indeed, for each 
trace class operator T, the map 6 ^ ti[TP Qe (X)} = tr[Te i0N P Q {X)e~ i0N ] is continuous)! 
and \ti[TE^ e (X)] \ < ||T||i||P^ e (X)|| < ||T||i, so the integral exists in the a-weak sense, and 
represents a bounded operator, with ||U fc (X)|| < 2n. 

Next, notice that for each k G N, the map X i— > V k (X) is an operator measure, that is, 
(T-additive in the weak operator topology. In fact, if (X n ) ne ^ is a sequence of mutually disjoint 
sets in £>(M), then for any I G N, 



i 

< Y,{v\P Qe {Xn)v) < (<p\p Qs (u™ =0 xm < y\ 



n=0 



-£e- ik9 (<p\P Qs (X n )<p) 

n=0 

so the dominated convergence theorem can be applied to get 

oo oo „ 2 7T r2n 

J2^\v k (xM = / e- ike (<p\pQ°(XM d6 = / e- ike (<p\pQ°(u™ =0 xM d6 

= ( V \V k (UZ X n ) V ). 

Let g : R — > M. be any bounded Borel function. Then the operator integral V k [g] = J g dV k 
can be defined in the cx-weak sense, as a bounded operator. This follows from the fact that g 
is bounded and 

|V^|(R) < 4 sup \tr[TV k (X)}\ < 87r||T||i 

for any trace class operator T, with \V^\ denoting the total variation of the complex measure 
= ti[TV k (-)]. (The map X \— > tr[TV rfc (X)] is a complex measure, because the weak and 
a-weak operator topologies coincide in a norm-bounded set.) 



3 This can easily be seen by e.g. considering a positive trace class operator T, using its spectral resolution 
and applying the strong continuity of the map 9 h- > e l6N . 



Lemma 5. For any bounded function g : R — > R, anc? a irace e/ass operator T, we have 

oo 

tr[TV*[</]] = 2tt ^(n|T|n + k) (n + k\g{Q)\n). 



n=0 



Proof. First note that (m\V k [g]\n) = tr[\n) (m\V k [g]] = J gdVJK,, by definition. On the other 
hand, for each X G B(R), we get 

U£ (X) = e- ike (m\P Q °(X)\n)d6= (m\P Q (X)\n) \ * e~ m e^'^dB . 



\n)(m\ ' 

</0 

This equals 27r(n + fc|P Q (X)|n) if m = n + and zero otherwise. Hence, 

(n + k\V k \g]\n) = 2v J g dP$ ){n+kl = 2n{n + k\g(Q)\n), 
and (m|y fc [^]|n) = whenever m ^ n + k. Thus, 

oo oo oo oo 

tr[TV k [g}] = ^(n\TV k [g]\n) = ^ ^(n\T\m)(m\V k [g}\n) = 2n^(n\T\n+k)(n+k\g{Q)\n) . 

n=0 n=0 m=0 n=0 

□ 

Now we are ready to prove the actual result of the paper. 

Theorem 1. Let S be a dense subset of [0, 2tt). The set {P® e \ 9 £ S} of observables is 
informationally complete. 

Proof. Let T G L(TC) be a selfadjoint trace class operator, such that tr[TP® e (X)} = for all 
X G <8(R) and 9 <E S. Since i— > tr[TP^"(X)] is continuous and S is dense it follows that 
t T [TP Q o(X)} = for all X G B(R) and G [0, 2tt). 

Let G N be fixed. By the definition of the operator measure V k , the assumption implies that 
V£{X) = tr[TV k {X)} = e- ike tr[TP Q e(X)]d9 = for all X G B(R). From the definition of 
the o"-weak integral V k [g] = f g dV k , it follows that tr[T^ fc [g]] = f gdV^ = for any bounded 
Borel function g : R — > R. Hence, by using Lemma [5l we get 

(4) ^2(n\T\n + k) (n + k\g{Q)\n) =0 

n=0 

for any bounded Borel function g : R — > R. We show by induction with respect to n that 
(n\T\n + k) = for all n G N. First, put g = /( fc ) in (jlj) (recall that / was a function satisfying 
the assumptions of Lemma [2]). By Lemma [U this gives (0|T|fc) = 0, which proves the initial 
step. The induction assumption is that for some m G N, (n\T\n + k) = for all n G N, n < m. 
We have to show that this implies (m + l|T|(ra + 1) + k) — 0. By the induction assumption 
and (pj, we have 

oo 

(n\T\n + k)(n + k\f( k+2 ( m+1 »{Q)\n) = 0, 

n=m+l 

where we have simply put g = /( fe+2 ( m+1 )), which is again a bounded Borel function. But, 
according to Lemma [H 

(n + k\f k+2(m+1 »{Q)\n) = 0, n>m + l, 
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and 

(( m + i) + k\f (k+2(m+1 »(Q)\m + 1)^0, 

so that necessarily (m + l\T\(m + 1) + k) — 0. This completes the induction proof. 

We have thus established that (l\T\l + k) = for all I, k G N. Since T is selfadjoint, this 
implies that T = 0, and the proof is complete. □ 

Note that the set S in the previous theorem can be chosen to be countable (e.g. S = 
[0, 2tt) HQ). Hence, in principle, it suffices to measure a sequence of quadrature observables in 
order to determine the state of the system. 

Though obvious, it may be worth to note that the quadrature observables Qg are not in- 
formationally complete in the sense of statistical expectation, that is, the numbers trfTQg], 
9 G [0,27r), do not, in general, determine the state T; for instance, the number states |n) are 
indistinguishable by the expectations, (n\Qg\n) = for all n and for all 9. 

4. WlGNER FUNCTION VS. PHASE SPACE DISTRIBUTIONS 

According to the result in the preceding section, the quadrature observables Qg, 9 G S (S C 
[0, 27r) is dense) constitute an informationally complete set of observables, i.e. the measurement 
statistics p® e , 9 G S, determine uniquely the state T of the quantum system. The question of 
experimental implementation of these observables is thus of utmost importance. 

The balanced homodyne detection with a strong auxiliary field is a well developed method of 
experimental quantum physics, and this method is known to yield the measurement statistics 
of the quadrature observable Qg, depending on the phase e ld of the (one-mode) auxiliary field. 
The heuristic physical argument behind this method is equally well known, see, e.g. [TQl [7] , the 
detailed mathematical justification being, however, more involved. 

If \z), z = re %6 , is the coherent state of the (one-mode) auxiliary field, then the actually 
measured observable in the balanced homodyne detection is given by a semispectral measure 
E z , whose first moment operator E z [l] is an extension of the restriction of the quadrature 
operator Qg on the domain D(a) of the signal mode operator a, and whose noise operator 
E z [2] — E z [l} 2 equals with the operator |r _2 iV (where N = a*a). Clearly, these results suggest 
that the high amplitude limit of E z is the spectral measure P Qb of Qg, notably since the spectral 
measures are known to be exactly those semispectral measures whose noise operators equal zero. 
There is, indeed, a rigorous quantum mechanical proof of the fact that in the high amplitude 
limit the observable E z tends to the spectral measure of Qg, though the actual meaning of this 
limit requires more caution [8j. 

The balanced homodyne detection scheme thus allows one to collect the statistics of the quad- 
rature observables. Assuming that the inverse Radon transform can be applied to the collection 
of distributions p^ e , 9 G [0, 27r), one obtains the Wigner function f^f of the state T, which is 
a unique representation of T. Clearly, there is no need to use the inverse Radon transform 
nor the Wigner function, since the set of quadrature observables is, in itself, informationally 
complete. In any case, this method of state reconstruction is unnecessarily complicated, since 
it requires that one measures infinite number of different observables Qg, 9 G S. 

It is well-known that a beam splitter combined with the Qo and Q n /2 -sensitive detectors 
constitutes a measurement of a phase space observable E D of the signal mode, with the gen- 
erating density operator D depending on the state of the (one-mode) auxiliary field, see, e.g. 
[U VII. 3. 7]. Using the balanced homodyne detection realization for the Q and Q w /2 -sensitive 
detectors one indeed obtains, in a rigorous quantum mechanical sense, an eight-port homodyne 
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detector realization of an arbitrary phase space observable E D [9]. A phase space observable 
E D is known to be informationally complete whenever the generating density operator D is 
such that tr[W gp D] 7^ for almost all phase space points (q,p) G M 2 ; here W gp is the Weyl 
operator, see, e.g. [2]. In particular, if the auxiliary field is idle, then the measurement statistic 
obtained from the eight-port homodyne detector is simply the Husimi function of the signal 
state T. Like the Wigner function, the Husimi function separates states, that is, the phase 
space observable E D , D = |0)(0|, is informationally complete. It is obvious that the state 
reconstruction via the eight-port homodyne detector is a huge simplification when compared 
with the reconstruction using only the balanced homodyne detection technique. 

Appendix: The proof of Lemma CD 

To begin the proof of the lemma, we first note that the Dawson's integral is clearly a C°°- 
function. We prove (a) by using the expansion 

, . 1 ^ (2?'-l)!! (2n-l)!' , > 

daw (x) = — > v ; , J + - —RJx), 

y J 2x^ (2x 2 y 2 n v ' 

where 



_ .2 OO OA 



x 2n ~ x f^jK 2 J -2n + l) 

and n = 1,2,3,... (see [HI equation 42:6:5, p. 407]). Since each derivative is either even or 
odd, it suffices to consider the limit x — > 00. Clearly, any derivative of the first part tends to 
zero at this limit, for any choice of n. As for R n , it is easy to see that for any given k, we get 
lim^oo R$h' (x) = for suffiently large n. Indeed, let 



Sn(y) = ^2~ 



y J 



o j\(2j-2n + iy 



so that Rn(x) = ^ 2n _i S n (x 2 ). Since the power series S n (y) is clearly convergent for all y G 
it can be differentiated k times (for any k G N) to get 

si k) (y) ^ 



^ (j-k)\(2j-2n + l) 

J—K 



From this we see that \Sn\y)\ < e y for any y > and n£N. Now for a fixed fceN, put e.g. 

n — k + 1. Since Rn\x) is clearly a finite sum of terms of the form -§z=r=i Sn \x 2 ), with 

—k < I < k, < I' < k, and A^y a constant, it follows that lim^^oo Rn\x) = 0. The proof of 
(a) is complete. 



To prove (b), consider the series 



1 (-l) n n! 



y ' 2^ 2 n (2n) v ' 

n=0 v ' 



We first note that the well-known relation 4-Hi(x) = 2lHi^i(x), I = 1, 2, 3 . . ., implies 

d m 

(6) — - H 2n {x) = 2n(2n-l)---(2n-m + l)2 m H 2n „ m (x), m < 2n. 

ax m 



Using the estimate \H n {x)\ < 'K2^ n y/nJ., where K > is a constant (P 22.14.17, p. 787]), 
we get 



'-l) n n\ d m TT , * 



2 n (2n)! cte r > 



< ^ K n\ 2n(2n - 1) • • • (2w - m + 1 V(2n - m)! ^ < gf g ggnl ^ 

(2n)! " ^(2ny. 



for all n, m G N, m < 2n. Now Y^=o T/^HiF < 00 ^ or an ^ m ^ ^ by the ratio test, which shows 

that the series obtained by differentiating ((5]) m times term by term converges uniformly in 
bounded intervals. Consequently, that series represents the mth derivative of the original series 
((51), the latter converging to some C^-function F : R — > R uniformly in bounded intervals. 
By using again the formula 4-Hi(x) = 2lHi-.\(x), we get 

n!(-l) n fT 2(n _ m) (x)_ 



(7) 2F^ m \x) = 2 2m J2^ 



n—m 

oo 



2 n (2(n-m))! ' 



(o\ 9F (2m+l)/ \ _ 9 2m+l n! ( l)"#2(n-m)-l 0*0 

W W ~ ^ 2»(2(n-m)-l)! " 

n=m+l v v y y 

In order to calculate the MacLaurin series of F, we need the expansion 
(9) (l-x)-(™ +1 ) = fV m + n V, -KKl, 

n=0 ^ ' 

which can be obtained from the binomial series (see 3.6.9. in [U p. 15]). 

Now, using @, the identity # 2 (n-m)(0) = (-l) n ~ m (2(n - m))!/(n - m)!, n > m (22.4.8 in 
P p. 777]), and © with x = |, we get 

2F (2 m)(0) = 2 2m (-l) m V — = 2 m (-l) m V ^— ^ 

n=m v ' n=0 



l) m m! ( n + m ) 4 = 2 2m+1 (-l) m m!. 

n=0 ^ ' 



Since F( 2m+1 '(0) = by © and the fact that if { (0) = for odd Z, we get the following 
MacLaurin series for F: 

do) f;(z£^ (2x) - I6H . 

m=0 V ' ' 

This is exactly the MacLaurin series for the first derivative of the Dawson's integral, since 



^ (-l) m m!2 2 
daw (^) = L (2m 



_2.2m+l 



m=0 



(see e.g. [H p. 406]). 

In order to prove that the series (fTUl ) actually converges to F (pointwise for all x G R) , we have 
to show that for each x G R, the corresponding remainder F^ k \^ k )x k / k\ , where G [— \x\], 
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goes to zero as k — > oo. By applying the estimate \Hi(y)\ < Ke^ y2 2^ 2 y/l\ to (CJ and (jHJ), we 
get 

2\F^(y)\ < e^2-Ky^±^ = Ke^2-m\y( n + m \^= 



2n)\ 

2\F^\y)\ < e ^^2»irf = ^e^2»( m + 1)! V (" + m + 

for all It is easy to see by induction that n!/ (2n) ! < n/2 n ~ l for all n > 1; using this, 

as well as the relation (m + l)2 m+1 = Xl^Lo ("^n™):^" which is obtained by differentiating (jHJ) 
and putting x = §, we conclude that 2\F^ m \y)\ < AKe^ y2 2 2m (m + 1)!, and 2\F i - 2m+l \y)\ < 
4V2Kei y2 2 2m+1 (m + 2)! for all y G R. Consequently, 

\F^(Qx 2m \ < 2 ^(2Nf>+l)! 
(2m)! ~ (2m)! 
| F(2m+ i) (0x2m+ i| ^ (2ixi) 2 -+ 1 (m + 2)! 

(2m + 1)! ^ V2Ae (2m + 1)! ' 

whenever i6l, ( G [— \x\] and m G N. It is easy to see that for each fixed x G K, the right 
hand sides of the inequalities tend to zero as m — > oo. Hence, we have shown that 

= V 1 J JT" (2x) 2m = daw (1) (x), x G R. 
z — ' (2ml! 

n=0 v ' 

This completes the proof of Lemma [U 
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